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Abstract 

We study a new approach to the problem of transparent boundary 
conditions for the Helmholtz equation in unbounded domains. Our ap- 
proach is based on the minimization of an integral functional arising from 
a volume integral formulation of the radiation condition. The index of 
refraction does not need to be constant at infinity and may have some 
angular dependency as well as perturbations. We prove analytical results 
on the convergence of the approximate solution. Numerical examples for 
different shapes of the artificial boundary and for non-constant indexes of 
refraction will be presented. 



1 Introduction 

Many problems in physical applications are modeled by wave propagation in un- 
bounded domains. Computing the numerical solution is a challenging problem 
and, usually, one has to introduce an artificial boundary to make the computa- 
tional domain finite. The problem of determining the boundary conditions on 
the artificial boundary has attracted many scientists and many methods have 
been studied leading to a very exciting field of mathematical research. 

In this paper, we consider time-harmonic electromagnetic wave propagation 
in unbounded domains. Let D C M. d , d > 2, be a (possibly empty) bounded 
domain. We consider the Helmholtz equation 

Au + k 2 n(x) 2 u = /, inR d \D, (1) 

where k > is called the wavenumber, n > is the index of refraction (which 
satisfies some assumptions to be specified later), / is a source term and u is part 
of the electromagnetic field. If D is not empty, boundary conditions on dD are 
given. To ensure the uniqueness of ([I]), it is necessary to introduce a radiation 
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condition at infinity. If n = 1 in the exterior of some ball, we can consider the 
following Sommcrfeld radiation condition 




— iku 




(2) 



where r = \x\ is the distance from the origin and i is the imaginary unit. 
The radiation condition ([2| follows from the expression of the far-field of the 
radiating solution. We will often refer as the solution of ([l])-([2]) as the outgoing 
solution. We shall discuss later the case when the index of refraction n is not 
constant outside any bounded domain. 

Depending on the assumptions on n, an exact formula for the solution of ([I]) 
may or may not exist (see for instance [il3"],|26]): however, when it exists, it may 
be intractable from a computational point of view. Thus, it is usually convenient 
to introduce an artificial computational bounded domain Q, with T = <9f2, and 
then impose some boundary conditions on T in such a way that the solution 

of the new problem approximates u in a good fashion, i.e. one has to try 
to eliminate spurious reflections of waves arising from the introduction of the 
artificial boundary T. 

As it is well-known, a large amount of work has been done on this kind 
of problems. The most used methods are based on local or nonlocal condi- 
tions involving u ([8].|14j.|19j). approximations of the Dirichlet-to- Neumann 
map ( [55] , [TS] , [5U] ) , infinite element schemes ([5J,[B]), boundary element meth- 
ods (see |33j ) and the perfectly matched layer method ( 7 ,[35 ). Many papers 
have been written of this subjects and for a deeper understanding of these 
problems and more recent developments, the interested reader can refer to 
[2], [6], [16], QI], P2], E2U, [22], [23], [34] and references therein. 

In this paper, we propose a new approach to the study of wave propagation 
problems in unbounded domains. Let us start from the simplest case of finding 
the outgoing solution u of the Helmholtz equation in M. d with constant index of 
refraction, i.e. to find u satisfying 



and the radiation condition ([2]). 

It is well known that other radiation conditions than ([2]) can be imposed at 
infinity to have the well-posedness of the problem. One of the most used is the 
following Rcllich- Sommcrfeld radiation condition 



which was proposed by Rellich in [3Qi . In the same paper, Rellich proved that a 
radiation condition can also be given in the form 



Condition Q can be considered as the starting point for our work, as we are 
going to explain shortly. 



Au + k 2 u = /, in R' 



d 



(3) 





(4) 
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A naive way to look at problem ([3])-(|4]) is the following: 
du 



Minimize 



Or 



iku 



dx among the solutions of Au + k 2 u = f in R . (5) 



Clearly, the outgoing solution is the unique minimizer of (§: it is the only 
solution of ([3]) that satisfies Q, since the integral in Q does not converge for 
any other solution of ph. 

Hence, when we introduce an artificial boundary 57, it is natural to consider the 
following constrained optimization problem: 



inf 



dv 
dr 



ikv 



1 

dx : Aw + k 2 v = f in ft >. 



(6) 



As the domain fl enlarges, it is reasonable to think that the minimizer of (§ 
approximates the solution of (|3])-(|4|- However, for a constant index of refraction, 
available methods in literature are very efficient and well-established. For this 
reason, in this paper we shall consider more general indexes of refraction. 

We now explain in details our method and the outline of the paper. We 
consider the Hclmholtz equation 



Au + k 2 n(x) 2 u = /, 



(7) 



here the index of refraction n is not necessarily constant at infinity but can have 
an angular dependency like 

n(x) — > noo{x/\x\) as \x\ —t +00. 

Under some appropriate assumptions on this convergence and on rioo, in [29] 
(see Theorem 1.4) the authors proved that a Sommerfeld condition at infinity 
still holds true under the form 



Wu{x) — ikn{x)u(x) - 



dx 



1+ x 



< +00. 



(8) 



There are other conditions of this type available in literature; we shall not discuss 
about them here and refer to , [35] , and references therein. 

Let O be a bounded domain and consider the following minimization problem 



inf{Jf2[u] : v satifies Av + k 2 n(x) 2 v = / in 57}, 



where 



Ja[v] 



\7v(x) — ikn(x)v(x)- 



dx 



(9) 



(10) 



The minimizer itn of Q is the approximation for the solution u of ([7])-(|8| that 
we propose in this work. 

Problem ^ can be easily implemented numerically: we have to minimize 
a quadratic functional in a bounded domain subject to a linear constrain. We 
shall prove that when the computational domain enlarges and tends to cover 
the whole space, the approximating solution un converges in H} oc norm to the 
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outgoing solution. However, when f2 is a ball and n is constant, most of the 
numerical methods cited before are more suitable and efficient for computing 
the solution than our approach. 

When the artificial domain O can not be chosen as a ball or as other domains 
with simple geometry, most of the methods available in literature may present 
some difficulty in their applications. Since we do not use any technique of sepa- 
ration of variables or analytic continuation, the implementation of our approach 
do not suffer the change of geometry of the artificial domain and we will show 
by numerical examples that such a change does not affect the computations. 

Our method can be applied also when the index of refraction is not constant 
outside any compact domain. In particular, we are thinking to an index of 
refraction that has some angularly dependency and may have some perturbation 
extending towards infinity. This is another relevant case in which our method 
is still applicable and of easy implementation. On the contrary, many of the 
methods cited before can not be applied (at least in a standard way). 

The paper is organized as follows. In Section [2] we prove a result of conver- 
gence for uq and other analytical results when the artificial domain is a ball. 
In Section [3] we give some numerical results on the convergence of the solution 
for a constant index of refraction for several choices of the artificial domain. 
In Section [4] we present some numerical examples for non-constant indexes of 
refraction. 



2 Mathematical framework and analytical re- 
sults 

In this section we assume that the artificial domain is a ball and prove results 
of convergence as the radius of the ball goes to infinity under quite general 
assumptions on n. To simplify the exposition, firstly we study equation ([I]) with 
D = and then the convergence result is generalized to the case D ^ with 
Dirichlet or Neumann boundary conditions on dD. 

Let U : M. d -> C be the solution of problem 0-([8]). In what follows, we shall 
assume that 

n£L°°(R d ) and n(x) > > for every x £ R d , (11) 

for some positive n*. We mention that the existence and uniqueness of the 
solution U is a challenging problem and they are usually obtained by using a 
limiting absorption principle. In this paper we will assume that U exists and 
is unique and refer to the papers cited in the Introduction for the appropriate 
assumptions on n. 

Let Br C M. d be the ball of radius R centered at the origin. For a solution u 

of 

Au + k 2 n(xfu = f in B R , (12) 



we mean a function u £ H 1 (B]i) satisfying the variational formulation of (12) 



j Vu-VTpdx + k 2 J n{x) 2 mpdx = J fpdx, \/ip£Cl{B R ). (13) 



Br Br 
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We denote by T> R the set of solutions of (13) in Bf 



V R = {uG H 1 {B R ) : u satifies (pi) in B R }. 



(14) 



Clearly, V R is a convex set. For u £ H 1 (B R ) we define the following seminorm: 

(15) 



Mr = 



Vu — iknu-, 



2 dx 
l + bl 



We consider the following family of minimization problems 

inf{[u] R : u € V R }, 

set 

d R = mf{[u] R : iieffl} 



(16) 



(17) 



and denote by ub r the minimizer of (161 (we will prove that ub r exists and is 
unique) . 



Lemma 2.1 Let u be a solution of 

Au + k 2 n(x) 2 u — in B R . 

Then 



Im / u r uda = 0. 



(18) 
(19) 



for every < p < R. 



Proof. This is a well-known lemma. It can be proved by multiplying (181 by 
"Xb p j integrating by parts and taking the imaginary part of the equation. □ 

The following lemma is crucial to prove H 1 convergence of the approximating 
solution. 



Lemma 2.2 Let u be a solution of (181 and assume that 

n n» ii 

w/iere 



njj(r) 



1 



I as 



uda, 



r\ JOB 



for every < r < R. Then we have 

(If 

(\Vu\' + k-n*\uy) — 

I Br, 



(\Vu\ 2 + k 2 n 2 \u\ 2 ) Y 
Proof. We assert that 



< 1- 1- 



21m / knu r u- 
b r 



dx 
'l + \x\ 



< 1 



(20) 
(21) 

(22) 
(23) 
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Indeed, we write 



Im 



dx 



knu r u- 

B R 1 + \X\ 



Im 



dx 



Im 



kn^u r u 



dx 



1 



\x\ 
(24) 



k(n — n»)u r u- — + 

1 + 1*1 

and notice that, since nj is a real- valued radial function, coarea formula and 

^(imf u r udv) 



qi9J> yield 
Im 



rf.r 



u r ud<r j dp = 0. 



(25) 



Since 



2fc(n — 7ij)M r u 



1+ d 



< 1 



(]u r \ 2 + k 2 n{x) 2 \u\ 2 )- 



dx 



+ \x\' 
(26) 



from (24)- (26 1 we easily get (23). From (23) and being 
K= J (|V U | 2 + fcVH 2 ) T ^ R -2Im 



knu r 



dx 
''l + \x\ 



we obtain ( 22 ) . 



□ 



In the following proposition, we prove that the minimizer of ( 16 1 exists, is 
unique and it is a solution of the Helmholtz equation in B R . 



Proposition 2.3 Let R > be fixed and assume that n satisfies ( 11 ) and ( 20 1 . 

Then, there exists a unique ub r G 2?r such that [ub r ]r — Ar, where d^ is given 

by m. 



Proof. Let {w m } m gN C be a minimizing sequence for (16) such that 

[u m ] R <d 2 R + ^, 
for any m £ N. The parallelogram identity yields 



(27) 



'Ui - u m - 


2 


'Ui + U m ~ 


2 


l 2 . 


+ 

R 


2 





[ut]% + [u m ] R 



for every <€N. Since (ue + u m )/2 € T>r : then 



< 



[ue} 2 R + [u m ] 2 R 



4- 



and from (27) we obtain 



[m - u m ] 2 R < 



for any I > m. From (22 ) and ( 11 ), we have that {u m } m gN is a Cauchy sequence 
in H 1 (B R ) and hence there exists u Rr & H (Br) such that u m — > ub r in 
H 1 (B R ). From such strong convergence we get that u Rr satisfies (13) and then 
u Br G V r . 

Now we prove the uniqueness by contradiction. Let us assume that there are 
two minimizers u and v. The parallelogram identity yields [u — v] R < which 



is a contradiction on account of Lemma 12.2 



□ 
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Lemma 2.4 Let U be the solution of ^ satisfying ^ and set 



Vw — iknu-, 



2 dx 

TTR 



The sequence {c?r}{.r>o} * s nondecreasing and d R < [U]oo, for every R > 0. 
Proof. Let R be fixed and consider p > 0. Since Ub r+p \ b G 2?r, we have 

hence, <i/{ is non-decreasing in R. Since U is a solution of the Helmholtz equation 
in R d , then C/| G V R , and thus d fl < [U) R < [U]^, for any R > 0. □ 

We are ready to prove our main theorem. 



Theorem 2.5 Let Rq > Q be fixed and assume that n satisfies (11) and (20 1 

Then 

lim ||u Bj? - C/|| ff i (B „)=0. 
Proof. Let p > R > Rq; since ub L G ^i? > w ^ have 



["Br - U Bp - 


2 

+ 


'u Br + u Bp - 


2 


i 2 


R 


2 


R 



Wb r }r + [ttB„]% 

2 



Since 2?r is a convex set, 



ub r + u Bp 



> df 



and, being [u Bp ] R < [u Bp ] p = d p , we have 



Wb r - u B ] R < 2{dl - 4). 



Lemma 2.2 and (111 yield 



/ (|V(u Bjl - u Bp )f + k 2 nl\u BR - u B f) -^r- < 2{d\ d R ), 

J B R 1 + Fl 



and, since Rq < R, 



\u Br -u Bp \\ H i iB ) < 



(d 2 p - d R ). 



We notice that Lemma |2.4| implies that d R converges to a finite number as 
R — » +oo and hence {u Br } r>Ro is a Cauchy sequence in H 1 (B Ro ). Thus Mb h 
converges strongly in H 1 (B Ro ) and the limit is a solution of (13 1 in B Ro . Since 
i?o is arbitrary, we can define a solution u of ( 13 ) in E d which satisfies (l8|). From 



the uniqueness of the solution, we get that u — U and the proof is complete. □ 



In the rest of this section, we show how to generalize Theorem |2.5| to problems 
in exterior domains. The generalization is straightforward but it needs some 
slight modifications as we are going to show in the following. 
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Let D c R d be a bounded piecewise regular domain containing the origin 
and consider the following problems: 



Am + k 2 n(x) 2 u = f, 
ju + (1 - j)u v = g, 

Vu — iknu- : — 



\D, 



on dD, 



2 dx 
1+la 



(28) 



< +oo, 



for j = 0,l,g e L 2 (dD) and where ^ is the outward normal to dD. We denote by 
U° D the solution of problem ( 28 ) (as before, we shall not specify the assumptions 



on n to get existence and uniqueness of the solution). 
For R > diam D, we consider the problems 



inl'^ / |Vu — iknux/\x\\ 2 ; — r : 

Ib r \d^ 1 + Fl 



Am + k 2 n(x) 2 u = /, in B R \ D, ju + (1 - j)u v = g, on dD }, (29) 



j = 0, 1. Let nj.c be given by 

% n D c \ Jo 







nda, if 9B r n L> c ^ 0, 
dB r nD" (30) 

if aB r n L> c = 0, 



r > 0, where D c =M. d \ D. We have the following result. 



Theorem 2.6 Let R > diam D be fixed and assume that n satisfies (111 and 
(20), where n* is replaced by ti^jj. Then there exists a unique minimizer u J D R 



of (29). Furthermore u 3 D R satisfies Au + k 2 n(x) 2 u = f in Br\D and we have 



D,R 

lim 1 1 u 

_R->+oo 



D, R _ U b\\Hi(B Ro \D) - 0) 



for any Rq > diam D. 

Proof. The proof follows the lines of the one of Theorem |2.5| However, since we 
are considering an exterior problem, Lemmas |2.1||2.2| must be adapted to this 
case. 

Let u be a solution of 

j Au + k 2 n(x) 2 u = 0, in Br \ £), 
}ju + (1 — j)u v = 0, on 3D. 

Let i?* be the radius of the largest ball contained in D centered at the origin. 
For R* < r < i?, we consider the domain 

r2 r = B r \ Z). 

If r > diam D, then it is easy to show that 

Im / u r uda = 0. 

JdB r 
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If i?* < r < diamD, then the assumptions on D guarantee that the divergence 
theorem holds on any connected component of Q r . From Helmholtz equation 
and the homogeneous boundary conditions on 3D, we obtain that 



Im 



r uda = 0, 



dB r nD" 



(31) 



for Rt, < r < R. Hence, we have proven the analogous of Lemma 2.1 for the 
exterior problem. To prove the analogous of Lemma |2.2| we have to take care 
of the application of the coarea formula in (25). This can be easily done by 



extending u and V« to zero in D. The rest of the argument needed to prove this 
theorem is completely analogous to what done before for the case D = 0. □ 



3 Numerical results for a constant index of re- 
fraction 

In this section, we study some convergence properties of the approximating 
solution given by the minimizer uq of {9j|. Our main goal in this section is to 
study the convergence for several choices of the shape of the artificial boundary 
in the case of constant index of refraction. 

We consider a classical model problem in scattering theory for which an 
exact solution is available. In particular, we consider 

'Au + k 2 u = Q, inM 2 \£^, 
u = g, on dBf , (32) 

Jr 2 \b- ~ ikux / \x\\ 2 dx < +oo, 

where Bf is the ball of radius f centered at the origin and g € L 2 (dB 7 ~) is a given 
function. We notice that, since n is constant, the radiation condition at infinity 



that we use in ( 32 ) is equivalent to the usual Sommerfeld radiation condition. 
This can be easily seen by expanding the solution at infinity in terms of negative 
powers of | x \ . We stress that in this case we use a functional without the weight 
(1 + | ^ej | ) 1 because we expect a better rate of convergence for the approximate 
solution. It is straightforward to verify that the results in Section [2] hold true 
also without such weight. 



As it is well-known, the exact solution to problem (32) can be obtained by 
separation of variables and is given by 

1 «-^°° f 27T Ff 1 (kr) ~ ~ ~ 

W(*,r) = -E fc0 / Q j^co s[ i(e -0 M 9)dO, 

where (9, r) denote the usual polar coordinates, H\ is the Hankel function of 
the first kind of order i and ^2 means that the term with £ = is multiplied by 
a factor 1/2. 

Let ft be a bounded domain containing Bf. For a function / 6 L 2 (<9f2) we 
consider the following problem: to find ua[f] such that 

'Au n [f} + k 2 u n [f] =0, mn\B f , 

PI'- <Mf}=9, on 8Bf 

C D u n [f] + C N d v u n [f] = f, on 50, 
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where Co and Cn may take the value or 1 according to the choice of the 
boundary condition to impose in dtt (Dirichlet or Neumann boundary condi- 
tion). The function / is the unknown of the problem and it must be chosen so 
that the corresponding solution uq[/] minimizes the integral 



AMf})= / Vu n [f](x)-iku Q [f}(x), 



dx. 



(33) 



n\B? 



We solve numerically the problem by using a (CGM) conjugate gradient 
method (see for instance |4j). In particular, we need to evaluate J (the Frechet 
derivative of J), and then consider the following two problems: 



!Ap + k 2 p = -f , mn\B f , ( AS[u] + k 2 5[u] = 0, in Cl\B f , 

p = 0, on dBf, P3 : i S[u] = 0, on 8B f , 

CnP + CNdvp = 0, on Oil, [Cz>S[u] + CNd„6[u] = w , on 

Problem P2 is the adjoint problem to Problem PI while Problem P3 com- 
putes the small variation S[u] of the solution during the iterative process of the 
CGM algorithm which is exposed below: 



Algorithm 1 Conjugate gradient algorithm 
Require: initial guess f = fo, tolerance e 

1: Solve PI 

2: Evaluate J (ttn[/]) 

3: Solve P2. 

4: Set go := p,wq := p, m = 1 

5: while stopping criteria 7 < e is not satisfied do 

6: Solve P3 

7-. Evaluate J (5[u]) 

8: Solve PI 

9: Set: p =< g m -i,w m -i > L i / < (C N p+ C D d u p),w m -i > L 2, 

g m = 9m-i - p(C N P + C D d v p) 
10: Set: u n [f] = u n [f] - pS[u], f = f - ptu TO -i 
11: Set: 7 =< g m ,w m -i > L 2 /||w m -i||| 2 
12: Set: w m = g m + jw m -i 
13: to := m + 1 
14: end while 



To solve Problems PI, P2, P3 we use a finite elements formulation; numerical 
simulations are performed through FREEFEM© software. The finest grid used 
contains up to 2 • 10 5 grid points, and we adapt the grid close to the artificial 
boundaries and where the large gradients of the solution form. Each simulation 
is firstly performed in a coarser grid and then the grid is gradually refined 
until a fully grid-independence is reached in the numerical results. We also have 
checked in which way the choice of the Dirichlet or Neumann condition on d£l 
affects the numerical results, and we have verified that this choice does not 
influence the final result in the case of smooth domains. However, Neumann 
boundary conditions are in general to prefer when one deals with coarse grid, 
as spurious oscillation are likely to manifest close to the boundary in the case 
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of Dirichlet boundary condition. For this reason, we present all the result by 
imposing Neumann condition on <9f2 (CV = 1, Cd = 0). 

We shall present numerical results for three choices of the artificial domain 
f2: a ball, an ellipse and a square. In each simulation e = 10~ 8 is the tolerance 
in Algorithm 1 and the Dirichlet boundary condition on dBf is given by g{6) — 
cos(j9), being j an integer. The radius of dBf is f — 1/2. 

The convergence properties of our method are tested through the evaluation 



of the usual error norms and the relative error of the integral (34), that is we 
compute: 



L 2 = (/ ||un[/]-<w ct || 2 dA) 1/2 , 

J A 

L 2 rel = L 2 /(f \\u exact \\ 2 dA)^ 2 , 

J A 

H 1 = ( f \\u n [f} - u exact \\ 2 + ||V(«n[/] - u exact )\\ 2 dA)^ 2 , 

J A 

Hid = H'/i f \\u exact \\ 2 + \\Vu exact )\\ 2 dA) 1 / 2 , 



A 



SJ re l = J(un[f]-U exact )/J(u exact ). 

We shall evaluate these errors in the whole domain fi\_Bx/2 and in the annular 
domain f2]y 2 = an d we check how the numerical solution converges to 

the exact one and how the shape of the artificial boundary of 51 affects this 
convergence. In the following paragraphs we shall present the results in the case 
in which fl is a ball, an ellipse, and a square. We have performed in each case 
the convergence analysis for the values j = 0, 2, 3, and k = 0.5, 1, 2. 



Balls The first case presented is relative to the annulus £lf/ 2 = Br\Bi/2- The 
numerical simulations are performed up to R = 8. In Tables [T] we report the 
errors on , 2 and fi^ 2 for R = 1, 2, 4, 8 (the Ll el ,Hl el errors on Q{/ 2 are also 
shown in Figs[l])a-f in lin-log scale). 

Ellipses The second case presented is relative to the domain E^ 2 = £r\Bi/ 2 
where £r is the ellipses of equation x 2 /(2R) 2 + y 2 /(R) 2 = 1. The numerical 
simulations are performed up to R — 8. We mention that, with respect to the 
previous case, the CGM method is slightly slower to reach the desired tolerance. 
In Tables[2]we report the errors on Q\^ 2 and E^ 2 for R = 1,2,4, 8 (the L 2 rel ,Hl el 

errors on il\/ 2 are a l so snown m Figsj2|a-f in lin-log scale). 

Squares We conclude this section by showing a numerical investigation on 
the domain S^ 2 = Qr\B 1 / 2 where Qr is the square centered in (0,0) and 
length side 2R. The numerical simulations are performed up to R — 8. Even if 
our approach can be applied to any domain, we mention that the presence of 
corners or edges deteriorate the convergence properties of the CGM method, at 
least if compared to the case in which the artificial boundary is regular (circle or 
ellipses) . We also notice that imposing the Dirichlet boundary condition on the 
problem PI (i.e. Cd = 1, Cm = 0) gives rise to a relevant oscillatory behavior of 
the solution near the corners due to the lack of regularity of the domain, while 
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Neumann boundary gives a more regular behavior in the numerical solution. 
In Tables |3j we report the errors on £l\/ 2 an d S±/2 °y varying the semi-length 
side R of Qr for R = 1, 2, 4, 8 (the L^. el ,H^ el errors on are also shown in 
Figs[3])a-f in lin-log scale). 



4 Numerical results for a variable index of re- 
fraction 

In this section we apply our approach to cases of variable indices of refraction, 
for which an exact solution is not generally given. We notice that most of the 
numerical methods available in literature can be of difficult application, since 
they usually need the knowledge of the exact solution in some exterior domain. 
We also mention that the perfectly matched layer might be of nontrivial im- 
plementation when the material is not an analytic function in the direction 
perpendicular to the boundary, see [24"]. [28]. 

We stress that for a variable index of refraction it can be hard to compute 
the asymptotic behavior of the integrand function in ^ and, unlike the case of 
constant index of refraction, the weight (1 + | | ) 1 in the integral (10) has to 
be maintained. 

Thus, given the same notation as the previous section, we are solving the 
following problem: to find / £ L 2 (dil) such that the solution ua[f] of 



Au n [f] + k 2 n(x) 2 un[f] = 0, 
un[f] = 9, on dB f 
d v u n [f]=f, on dn, 



in d\Bf 



minimizes the integral 



Vu n [f}(x) - ikn(x)ua[f]{x) 



n\B f . 



:dx. 



(34) 



In the following paragraphs we shall propose two indexes of refraction, 
namely 



ni(x,y) 
n%{x,y) 



2 + (exp(-(z - l) 2 - y 2 ) + exp(-(x + l) 2 - y 2 )) 



2 + a 



\/x 2 + y 2 ' 



\/x 2 + y 2 ' 



being a a non negative real number. These indexes are shown in Figs]4^-b, with 
a = 0.1. We stress that ri\ has an exponential decay to 2 at infinity and that n% 
it is not constant at infinity. For both cases, the geometry of the computational 
domain is the annulus Qy 2 = ^R\B\/2 for which a better convergence rate is 
achieved in the CGM algorithm, as seen in the previous section for a constant 
index of refraction. 

Being an exact solution not available, we test the convergence properties 
as follows: we set the solution for R = 16 as the reference one, and we shall 
compute the error norms up to R = 8 as done in the previous sections. We 
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Index of refraction n± (x) 



Index of refraction (a:) 





(a) ni (x, y) (b) nJj- 1 ^, j/) 

Figure 4: The indexes of refraction rix and n 2 . 
always consider j = 0, k = 1. 



We notice that the choice of the functional (34) is appropriate only when 



the index of refraction does not generate guided modes (trapped waves). Due 
to the exponential decay at infinity, this requirement is satisfied by ri\ (see [1 
and [SI]). Regarding the index n 2 it is known from that there is a critical 



value 



according to which no guided modes are generated for a < a 9 " 1 . It 



is not easy to give a quantitative estimate of a gm . By evaluating the integral 
( 34 1 for several values of a (see Fig{5] where the integral is shown at various 
exterior radius R in log-lin scale), it seems that for a ^ 0.6 the integral diverge 
as R — > oo. This suggests to consider a ;$ 0.6 for which it appears that the 
functional is bounded as R — > co. The numerical results for n 2 are therefore 
shown only for the small value a = 0.1: if we are luckily in the case in which 
no guided modes are generated, our approach works well and the numerical 
solution converges to the exact one. More accurate estimates on a gm will be the 
object of future work. 

Case n 1 (x): In Tables |4j we report the errors on and £ly 2 ^ or ^ = 

1, 2, 4, 8, (the L^. el ,H* el errors on fij , 2 are also shown in Fig|7|)a) in lin-log scale. 
In Figs 6^,-b the real and imaginary parts of the numerical solution Mjjk [/] are 
shown for R = 16 

Case ri2 A {x): In Tables jHJ we report the errors on fl\/ 2 an d f° r R = 

1, 2, 4, 8 (the L^. el ,H^ el errors on fi^ 2 are also shown in Figj9j)b) in lin-log scale. 
In Figs 8^i-b the real and imaginary parts of the numerical solution Uqk [/] are 
shown for R = 16. 
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Value of the integral J(u^r ), n(x) = n^ix) 




Figure 5: The integral value J(uqr yvj), by varying the exterior radius R for 
various a (log-lin scale): for a ^ 0.6 a likely numerical divergence is noticeable. 





-5 5 



-15 -10 



(a) 



(b) 



Figure 6: The real and imaginary parts of the numerical solution Uqw , index of 

1/2 

refraction ni(x,y). 
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-10 -5 5 10 15 -15 -10 -5 5 10 15 



(a) n2' 1 (x,y) (b) r&^ix, y) 

Figure 8: The real and imaginary parts of the numerical solution Mqi6 , index 

1/2 

of refraction n!] 1 ^,?/). 
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relative error on f3[/ a for index of refraction n"' 1 
-2.5 1 1 1 1 1 1 




1 2 3 4 5 6 7 8 



R 



Figure 9: Numerical errors (lin-log scale) of the solution for variable index of 
refraction n^ 1 on the domain £l\j 2 by varying the radious R of ^i/ 2 - 

Conclusions 

We studied a new approach to the problem of transparent boundary conditions 
for the Helmholtz equation. Our method is based on the minimization of an 
integral functional which arises from a volume integral formulation of the Som- 
merfeld radiation condition. When the artificial domain is a ball, we proved 
analytical results on the convergence of the approximate solution to the exact 
one under quite general assumptions on the index of refraction. We tested nu- 
merically our method for a constant index of refraction by choosing three shapes 
of the artificial boundary in R 2 : a ball, an ellipse and a square. The convergence 
results are satisfactory; in particular we observed a good convergence rate for 
the H^ el norm. We also applied our approach to other indices of refraction, one 
exponentially decaying at infinity to a constant, and one not constant at infinity. 

Since our approach is of easy implementation and has wide application, we 
believe that it can be a good starting point for studying problems with non- 
trivial indices of refraction. However, when the index of refraction is constant, 
the methods available in literature are more suitable. We stress that numeri- 
cal results show that our approach do not suffer the choice of the shape of the 
artificial domain. 

Our approach applies also to more general settings. In particular, we have in 
mind the waveguide problem. In this case, the functional J has to be modified 
on account of the presence of guided modes (see [2], [9] -[32]); this will be the 
object of future work. 



19 



References 

[1] S. Agmon, Spectral properties of Schrddinger operators and scattering 
theory. Ann. Sc. Norm. Sup. Pisa Ser., 4 2 (1975), 121-218. 

[2] O. Alcxandrov and G. Ciraolo. Wave propagation in a 3-D optical 
waveguide. Math. Models Methods Appl. Sci. (M3AS), 14 (2004), no.6, 
819-852. 

[3] R.J. Astley, Infinite elements for wave problems: a review of current 
formulations and an assessment of accuracy, Int. J. Numer. Methods 
Engrg., 49 (2000), 951-976. 

[4] C. Atamian, Q.V. Dinh, R. Glowinski, J. He, J. Pcriaux Control ap- 
proach to fictitious domain methods: Applications to fluid dynamics and 
electromagnetics. Proc. Fourth Int. Conf. on Domain Decomposition 
Methods for Partial Differential Equations, SIAM, Philadelphia (1991), 
pp. 275-309. 

[5] P. Bcttcss, Infinite elements, Int. J. Numer. Methods Engrg., 11 (1977), 
53-64. 

[6] P. Bettess, Infinite Elements, Pcnshaw Press, Sunderland, UK, 1992. 

[7] J. -P. Berenger, A perfectly matched layer for the absorption of electro- 
magnetic waves, J. Comput. Phys. 114 (1994), 185-200. 

[8] A. Bayliss and E. Turkel, Radiation boundary conditions for wave-like 
equations, Comm. Pure Appl. Math., 33 (1980), 707-725. 

[9] G. Ciraolo, A method of variation of boundaries for waveguide grating 
couplers, Applicable Analysis, 87 (2008), 1019-1040. 

[10] G. Ciraolo, A radiation condition for the 2-D Helmholtz equation in 
stratified media, Comm. Part. Diff. Eq., 34 (2009), 1592-1606. 

[11] G. Ciraolo and R. Magnanini. Analytical results for 2-D non-rectilinear 
waveguides based on the Green's function. Math. Methods Appl. Sci., 
31 (2008), no.13, 1587-1606. 

[12] G. Ciraolo and R. Magnanini, A radiation condition for uniqueness in 
a wave propagation problem for 2-D open waveguides, Math. Methods 
Appl. Sci., 32 (2009), 1183-1206. 

[13] D. Colton and R. Kress, Inverse Acoustic and Electromagnetic Scatter- 
ing Theory, Springer- Verlag Berlin, Heidelberg, New York, 1992. 

[14] B. Engquist and A. Majda, Radiation boundary conditions for acoustic 
and elestic wave calculations, Comm. Pure Appl. Math., 32 (1979), 314- 
358. 

[15] D. Givoli and J.B. Keller, A finite element method for large domains, 
Comput. Methods Appl. Mech. Engrg., 76 (1989), 41-66. 



20 



D. Givoli, Numerical Methods for Problems in Infinite Domains. Stud- 
ies in Applied Mechanics, vol. 33, Elsevier Scientific. Publishing Co., 
Amsterdam, 1992. 

D. Givoli, Exact representations on artificial interfaces and applications 
in mechanics, AMR, 52 (1999), 333-349. 

D. Givoli, Recent advances in the DtN FE method, Arch. Comput. Meth- 
ods Engrg., 6 (1999), 71-116. 

D. Givoli, High-order local non-reflecting boundary conditions: a review, 
Wave motion, 39 (2004), 319-326. 

M.J. Grote and J.B. Keller, On nonreflecting boundary conditions, J. 
Comput. Phys., 122 (1995), 231-243. 

T. Hagstrom, Radiation boundary conditions for the numerical simula- 
tion of waves, Acta Numer., 8 (1999), 47-106. 

I. Harari, A survey of finite element methods for time-harmonic acous- 
tics, Comput. Methods Appl. Mech. Engrg., 195 (2006), 1594-1607. 

F. Ihlenburg, Finite Element Analysis of Acoustic Scattering, Springer- 
Verlag, New York, 1998. 

S.G. Johnson, Notes on perfectly matched layers (PMLs). Massachusetts 
Institute of Technology, Technical Report, 2007; updated 2010. Avail- 
able at:http://www-math. mit.edu/stevenj/18.369/pml.pdf. 



J.B. Keller and D. Givoli, Exact nonreflecting boundary conditions, J. 
Comput. Phys., 82 (1989), 172-192. 

R. Leis, Initial Boundary Value Problems in Mathematical Physics, J. 
Wiley & Teubner Verlag, Stuttgart, 1986. 

M. Medvinsky, E. Turkel, U. Hetmaniuk, Local absorbing boundary con- 
ditions for elliptical shaped boundaries. J. Comput. Phys. 227 (2008), 
no. 18, 8254-8267. 

A. F. Oskooi, L. Zhang, Y. Avniel, S.G. Johnson, The failure of perfectly 
matched layers, and towards their redemption by adiabatic absorbers. 
Opt. Express, 16 (15) (2008), pp. 11376-11392. 

B. Pcrthame and L. Vega, Energy concentration and Sommerfeld condi- 
tion for Helmholtz equation with variable index at infinity, GAFA Geom. 
Funct. Anal., 17 (2008), 1685-1707. 

F. Rellich, Uber das asymptotische Verhalten del Losungen von Au + 
Xu — in unendlichen Gebieten, Jber. Deutsch. Math. Verein., 53 
(1943), 57 - 65. 

Y. Saito, Spectral and scattering theory for second-order differential op- 
erators with operator-valued coefficients, Osaka J. Math. 9 (1972), 463- 
498. 



21 



[32] Y. Saito, Schrodinger operators with a nonspherical radiation condition, 
Pacif. J. of Math., 126 (1987), 331-359. 

[33] S. A. Sauter and C. Schwab. Boundary element methods. Springer Series 
in Computational Mathematics, 39. Springer- Verlag, Berlin, 2011. 

[34] S.V. Tsynkov, Numerical solution of problems on unbounded domains. 
A review, Appl. Numer. Math., 27 (1998), 465-532. 

[35] E. Turkel and A. Yefet, Absorbing PML boundary layers for wave-like 
equations, Appl. Numer. Math. 27 (1998), 533-557. 

[36] M. Zubeldia, Energy concentration and explicit Sommerfeld radiation 
condition for the electromagnetic Helmholtz equation, preprint 2012, 
arXiv: 1201.0494. 



22 



Constant index of 
Errors 



. on Si} 



rcfr action: 
(-10- 1 ) 



Constant index of refraction: 
Errors o n Of^QlO" 1 ) 



el 



H 



H 



0.5 



0.664 
0.341 
0.161 
0.113 



0.509 
0.259 
0.123 
0.0865 



2.06 
1.05 
0.498 
0.352 



1.21 
0.615 
0.292 
0.205 



3.22 
1.21 
0.854 
0.499 



1.61 
1.95 
1.72 



0.679 

0.52 

0.308 



2.06 
2.07 
1.97 



0.714 
0.471 
0.309 



0.421 

0.211 

0.139 

0.0661 

0.261 

0.138 

0.0788 

0.0178 



0.331 

0.163 

0.107 

0.0512 

0.205 

0.108 

0.0621 

0.0141 



1.19 

0.642 

0.423 

0.201 

0.767 

0.407 

0.232 

0.0526 



0.602 

0.308 

0.203 

0.0967 

0.251 

0.133 

0.0759 

0.0172 



2.25 

1.17 

0.941 

0.452 

0.112 

0.836 

0.485 

0.112 



0.851 
0.801 
0.589 

0.318 
0.265 
0.0933 



0.374 
0.226 
0.112 

0.143 

0.0777 

0.0186 



1.09 

0.121 

0.841 

0.699 
0.596 
0.232 



0.315 
0.231 
0.111 

0.135 

0.0805 

0.0179 



1.88 
0.211 
0.0209 
0.00327 



3.41 
0.383 
0.0379 
0.00595 



7.77 
0.872 
0.0863 
0.00135 



3.13 
0.351 
0.0348 
0.00545 



0.582 
0.158 
0.0654 
0.00168 



1.68 
1.23 
0.693 



2.96 
1.85 
0.993 



3.71 
1.58 
0.753 



1.43 

0.609 

0.289 



1.86 
0.245 
0.0431 
0.0218 



0.323 
0.425 
0.0749 
0.0379 



7.62 
1.01 
0.176 
0.895 



3.09 
0.407 
0.0716 
0.0363 



0.687 
0.168 
0.0134 
0.107 



1.48 
1.21 
0.878 



2.69 
1.49 
0.905 



3.62 
1.59 
1.21 



1.38 

0.596 

0.435 



1.66 
0.357 
0.185 
0.102 



2.51 
0.539 
0.281 
0.155 



6.59 
1.41 
0.736 
0.406 



2.61 
0.559 
0.291 
0.161 



0.905 
0.0836 
0.0282 
0.0172 



1.59 
1.21 
0.924 



1.61 

0.933 

0.519 



3.13 
2.58 
2.12 



1.03 

0.715 

0.464 



1.56 
0.0446 
0.00113 
0.000051 



3.63 
0.103 
0.00263 
0.000118 



8.81 
0.232 
0.00591 
0.000266 



2.66 
0.0757 
0.00192 
0.0000865 



0.0143 
0.00951 
0.000453 
0.000021 



0.721 
0.252 
0.0957 



1.618 
0.564 
0.214 



1.88 

0.387 

0.108 



0.608 
0.124 
0.0348 



1.61 

0.0569 

0.00288 

0.000921 



20573 
0.131 
0.00658 
0.00211 



8.33 
0.294 
0.0149 
0.00477 



2.71 
0.0959 
0.004816 
0.00155 



0.0561 
0.0351 
0.000396 
0.000216 



0.822 
0.352 
0.211 



0.181 
0.763 
0.452 



0.208 
0.511 
0.299 



0.673 
0.164 
0.0963 



1.73 
0.121 
0.0415 
0.0181 



3.65 
0.255 
0.0878 
0.0383 



8.71 
0.609 
0.209 
0.0913 



2.83 
0.198 
0.0681 
0.0297 



0.367 
0.0654 
0.00143 
0.000956 



1.09 

0.707 

0.456 



2.04 
1.19 
0.661 



2.55 
1.58 
1.09 



0.808 
0.494 
0.305 



Table 1: The L 2 , Lf el , H^ el ,SJ re i errors of the numerical solution on the annulus 
f2]y 2 and ^f/ 2 by varying the exterior radius R of the annular domain fif/ 2 for 



1/2 

the case j = 0, 2, 3, k = 0.5, 1, 2. 
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Constant index of 
Errors on 



refraction: 

(•lo- 1 ) 



Constant index of 
i on E* /2 



Er 



refraction: 
I-IO" 1 ) 



el 



H 1 



Hi 



S.J r 



el 



H 1 



Hi 



0.5 



0.572 
0.291 
0.142 
0.0798 



0.435 
0.221 
0.108 
0.0607 



1.781 
0.902 
0.441 
0.245 



1.04 
0.528 
0.258 
0.144 



2.11 
0.966 
0.49 
0.581 



1.27 
2.24 
2.27 
1.54 



0.711 
0.761 
0.503 
0.232 



2.31 
2.646 
2.51 
1.731 



1.03 
0.751 
0.478 
0.229 



0.382 
0.179 
0.0905 
0.0496 



0.296 
0.139 
0.0701 
0.0384 



1.19 
0.553 
0.275 
0.151 



0.574 
0.265 
0.132 
0.0726 



2.78 
1.15 
0.641 
0.231 



0.0807 
1.061 
0.671 
0.549 



0.0467 
0.379 
0.159 
0.0881 



1.55 
1.45 
0.967 
0.76 



0.576 
0.345 
0.156 
0.0863 



0.5 



0.219 

0.115 

0.0585 

0.00601 

1.282 

0.161 

0.0239 

0.00964 



0.172 

0.0907 

0.0461 

0.00474 

2.232 

0.292 

0.0435 

0.0175 



0.689 

0.342 

0.173 

0.00178 

5.12 

0.548 

0.0795 

0.0305 



0.225 

0.112 

0.0523 

0.00582 

2.051 

0.221 

0.0321 

0.0123 



0.432 

0.211 

0.115 

0.0371 

0.449 

0.0541 

0.00249 

0.00125 



0.405 

0.332 

0.249 

0.045 

2.18 

1.52 

1.25 

0.715 



0.241 

0.122 

0.0611 

0.00772 

3.71 

2.35 

1.81 

0.991 



0.944 

0.749 

0.559 

0.101 

6.42 

2.37 

1.51 

0.794 



0.237 

0.119 

0.0603 

0.0076 

2.53 

0.921 

0.579 

0.303 



1.27 
0.0186 
0.0371 
0.0261 



2.21 
0.324 
0.0658 
0.0451 



5.05 
0.679 
0.133 
0.0921 



2.05 
0.275 
0.0541 
0.0372 



0.525 
0.072 
0.0211 
0.00925 



2.21 
1.97 
0.844 
0.842 



3.47 
2.61 
0.951 
0.775 



6.45 
3.31 
1.17 
1.09 



2.54 
1.25 
0.432 
0.32 



1.14 
0.325 
0.168 
0.00182 



1.73 
0.491 
0.253 
0.00275 



4.51 
1.24 
0.636 
0.00711 



1.78 
0.492 
0.251 
0.00279 



0.634 
0.0211 
0.0132 
0.0026 



2.04 
1.76 
1.19 
0.0291 



2.54 
1.58 
0.781 
0.0141 



6.22 
3.81 
2.66 
0.0661 



2.25 
1.16 
0.652 
0.0127 



0.5 



0.942 
0.0361 
0.00251 
0.000482 



2.18 
0.0841 
0.00581 
0.00112 



4.51 
0.127 
0.00849 
0.00163 



1.46 
0.0412 
0.00272 
0.00053 



0.0831 
0.00354 
0.00249 
0.00000647 



1.44 
0.394 
0.145 
0.0895 



3.28 
0.884 
0.326 
0.201 



5.83 
0.729 
0.178 
0.101 



1.83 
0.235 
0.0575 
0.0324 



0.887 

0.0584 

0.00471 

0.00265 

1.02 

1.04 

0.0415 

0.000811 



2.02 

0.133 

0.0107 

0.00605 

2.16 

0.221 

0.0878 

0.0171 



4.19 

0.217 

0.0168 

0.00921 

4.98 

0.461 

0.178 

0.0272 



1.36 

0.0707 

0.0054 

0.00299 

1.62 

0.151 

0.0581 

0.000881 



0.141 

0.0112 

0.00127 

0.000201 

0.321 

0.0201 

0.00201 

0.00000286 



1.42 

0.774 

0.281 

0.185 

1.82 

1.17 

0.638 

0.0178 



3.18 

1.68 

0.603 

0.392 

3.64 

2.09 

1.07 

0.0237 



5.57 

1.56 

0.394 

0.258 

7.03 

2.71 

1.41 

0.0401 



1.81 

0.503 

0.126 

0.0831 

2.25 

0.851 

0.431 

0.0118 



Table 2: The L 2 7 Ll el ,H^ el ,SJ re i errors of the numerical solution on f2^ 2 and 
Ey 2 by varying the semi-minor axis R of E^ 2 for the case j = 0, 2, 3, k = 



0.5,1,2. 
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E 



rrors on 
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Constant index o 
Errors on S^. 



refraction: 
(•10- 1 ) 



el 



H 1 



Hi 



8J r 



el 



H 1 



Hi. 



0.631 
0.321 
0.153 
0.101 



0.481 
0.243 
0.116 
0.0771 



1.95 
0.986 
0.471 
0.312 



1.14 
0.577 
0.276 
0.182 



2.81 
0.571 
0.231 
0.0721 



0.879 
1.84 
2.07 
1.65 



0.594 
0.717 
0.516 
0.278 



2.14 
2.269 
2.22 
1.88 



1.13 
0.727 
0.476 
0.277 



0.0421 
0.197 
0.129 
0.0467 



0.327 
0.153 
0.101 
0.0362 



1.302 
0.603 
3.93 
0.142 



0.625 
0.289 
0.181 
0.0681 



0.806 
0.606 
0.431 
0.0832 



0.581 
0.944 
0.701 
0.489 



0.402 
0.381 
0.201 
0.0881 



1.431 

1.23 

1.181 

0.684 



0.621 
0.331 
0.212 
0.0854 



0.5 



0.241 

0.139 

0.0497 

0.0320 

1.745 

0.193 

0.0192 

0.00315 



0.188 

0.109 

0.0391 

0.0251 

3.17 

0.351 

0.0351 

0.00572 



0.729 

0.411 

0.146 

0.03942 

7.26 

0.799 

0.0796 

0.0131 



0.238 

0.134 

0.0479 

0.0308 

2.92 

0.322 

0.0321 

0.00525 



1.935 

0.851 

0.321 

0.0851 

0.594 

0.143 

0.0675 

0.0016 



0.317 

0.336 

0.207 

0.176 

2.01 

1.91 

1.25 

0.695 



0.223 

0.141 

0.0561 

0.0332 

3.59 

3.02 

1.88 

0.993 



0.798 

0.798 

0.473 

0.395 

8.53 

3.96 

1.65 

0.764 



0.235 

0.143 

0.0571 

0.0322 

3.41 

1.53 

0.636 

0.293 



1.73 
0.227 
0.0415 
0.0141 



3.07 
0.394 
0.072 
0.0245 



7.14 
0.931 
0.171 
0.0578 



2.89 
0.377 
0.0691 
0.0234 



0.681 
0.149 
0.0133 
0.0079 



1.98 
1.94 
1.22 
0.654 



3.39 
2.74 
1.48 
0.663 



8.39 
3.91 
1.66 
0.904 



3.35 
1.491 
0.618 
0.323 



1.51 
0.344 
0.139 
0.0961 



2.37 
0.52 
0.211 
0.144 



6.25 
1.365 
0.553 
0.381 



2.47 
0.539 
0.218 
0.151 



0.912 
0.0624 
0.0144 
0.0852 



1.76 
1.62 
1.02 
0.891 



2.587 
1.67 
0.771 
0.492 



7.08 
3.37 
2.19 
2.01 



2.71 
1.10 
0.59 
0.435 



0.5 



1.17 
0.0422 
0.00237 
0.000191 



2.71 
0.0981 
0.00551 
0.000445 



6.13 
0.188 
0.00861 
0.000664 



1.99 
0.0611 
0.00279 
0.000216 



0.97 
0.0132 
0.00292 
0.000112 



1.63 
0.704 
0.247 
0.0962 



3.75 
1.58 
0.554 
0.215 



7.56 
1.69 
0.393 
0.109 



2.45 
0.546 
0.115 
0.0351 



1.21 

0.0504 

0.00304 

0.00108 

0.541 

0.0988 

0.0276 

0.0187 



2.76 

0.115 

0.00695 

0.00246 

1.14 

0.208 

0.0583 

0.0396 



6.31 

0.234 

0.0137 

0.00437 

2.73 

0.497 

0.138 

0.0923 



2.05 

0.0761 

0.00445 

0.00142 

0.891 

0.162 

0.0451 

0.0301 



0.148 

0.0201 

0.00986 

0.000146 

0.337 

0.0401 

0.000102 

0.000405 



1.69 

0.824 

0.356 

0.172 

0.929 

1.12 

0.501 

0.478 



3.81 

1.81 

0.771 

0.367 

1.91 

2.07 

0.931 

0.673 



7.76 

1.92 

0.521 

0.241 

3.42 

2.63 

1.25 

1.07 



2.51 

0.621 

0.167 

0.0777 

1.21 

0.833 

0.392 

0.323 



Table 3: The L 2 7 Ll el ,H^ el ,SJ re i errors of the numerical solution on f2^ 2 and 

Sy 2 by varying the semi-lenght side R of S^ 2 for the case j — 0, 2, 3, k = 
0.5,1,2. 



25 











Variable index of refraction ni 


(x): 




Variable index of refraction m(x): 














Errors 


on Q\ /2 


(•lo- 1 ) 








Errors on £~2^y 2 


(-10- 1 ) 






3 


k 


R 




L 2 




H 1 


Hl el 


SJ re , 




L 2 




H 1 


HL 


SJrel 







1 


1 




0.734 


0.582 


2.85 


1.01 


0.088 




















2 




0.284 


0.225 


1.12 


0.403 


0.156 




1.21 


0.559 


1.93 


0.499 


0.672 








4 




0.174 


0.137 


0.502 


0.177 


0.0511 




0.821 


0.242 


1.49 


0.229 


0.157 








8 




0.0656 


0.0521 


0.203 


0.0728 


0.00602 




0.411 


0.0826 


0.791 


0.0842 


0.00193 





Table 4: The L 2 , L 2 el , Hl ell 5J re i errors of the numerical solution on the annulus 
Q\/ 2 an d ^i/2 Dv var y m g the exterior radius R of the annular domain £lf/ 2 for 
the variable index of refraction ni(x) — 2 + (exp(— (x — l) 2 — y 2 ) + exp(— (x + 
l) 2 -y 2 ))cos(#). 











Variable index of refraction n^" 1 (x): 




Variable index of refraction n®' 1 (x): 












Errors on Q\ /2 (-10 -1 ) 








Errors 


on nf /2 (-10- 1 ) 




3 


k 


R 




L 2 




H 1 




5J rel 




L 2 


T 2 


H 1 




6J rB i j 





1 


1 




0.522 


0.408 


1.65 


0.582 


0.826 


















2 




0.221 


0.179 


0.715 


0.249 


0.442 




0.531 


0.237 


1.03 


0.238 


0.0925 






4 




0.121 


0.0951 


0.402 


0.141 


0.273 




0.421 


0.122 


0.0801 


0.128 


0.0833 






8 




0.0872 


0.0681 


0.245 


0.0861 


0.168 




0.478 


0.0951 


0.911 


0.101 


0.0161 



Table 5: The L 2 , L 2 ell H} ell 5J re i errors of the numerical solution on the annulus 
f2]y 2 and ^ v varvm g the exterior radius R of the annular domain Qy 2 f° r 
the variable index of refraction n^ix) = 2 + 0.1 cos(6*) 
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